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Q> , The one-hole spectral weight for two chains and two dimensional lattices is 



S 



X 



studied numerically using a new method of analysis of the spectral function 
within the Lanczos iteration scheme: the Lanczos spectra decoding method. 



Q , This technique is applied to the t — Jz model for — > 0, directly in the infinite 

O 

size lattice. By a careful investigation of the first 13 Lanczos steps and the 

first 26 ones for the two dimensional and the two chain cases respectively, we 

get several new features of the one- hole spectral weight. A sharp incoherent 

peak with a clear momentum dispersion is identified, together with a second 

broad peak at higher energy. The spectral weight is finite up to the Nagaoka 

energy where it vanishes in a non-analytic way. Thus the lowest energy of one 

hole in a quantum antiferromagnet is degenerate with the Nagaoka energy in 

the thermodynamic limit. 
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After the discovery of superconductivity in materials which exhibit antiferromagnetic 
(AF) long range order in the insulating phase, much attention has been given to the study 
of strongly correlated system at low doping. Despite the enormous amount of work, the 
dynamical properties of one hole in a quantum antiferromagnet (QAF) are still a subject of 
debate 

A satisfactory model of a QAF is provided by the Heisenberg hamiltonian which, in 2D, 
is known to show Neel long range order. In the following we mimic the properties of a QAF 
by use of the simpler AF Ising hamiltonian whose exact ground state is the Neel state and 
we address the problem of the propagation of a hole in the Jz —>■ limit of the t — Jz model 
defined by: 

H=-t y: (4c,. + h.c.) +Jzy: St s] (1) 

<i,j>,a <i,j> 

where the constraint of no double occupancy is understood. Both the two chain (2C) and 
the two dimensional square lattice (2D) will be considered. The Nagaoka (NK) theorem 

provides a lower bound ep to the ground state energy of one hole in the — > limit: 
ep = —zt, where z is the coordination number of the lattice {z = 3 for the 2C problem and 
z = 4 for the 2D case). 

A complete description of the one particle spectrum, not limited only to the ground 
state, was first given in the seminal work of Brinkman and Rice (BR), where the so called 
"retraceable path approximation" (RP) was introduced ||5[. In this approximation the spec- 
tral weight is completely incoherent and the hole is essentially localized without dispersion. 
The RP is exact in ID ( where NK theorem does not apply) and recently it has also been 
shown to be exact in the limit of infinite spatial dimensionality . However, no analytical 
solution is available for finite z > 2 and it is not clear how accurate the RP is in these cases. 
Indeed, in finite dimensions D > 1 several one hole paths exist, allowing for the propagation 
of the hole §. 

Recently, the problem of hole propagation has been addressed by use of the quite efficient 
Lanczos technique in 2D lattices with up to 26 sites [0 or with a promising approximate 
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technique up to 50 sites ||^. Despite the remarkable numerical effort, still the number of 
sites appears to be too small for a precise determination of the properties of the spectral 
function or for the ground state properties for — > 0. 

In this work, we present a numerical study of the one hole spectral weight by means of 
the Lanczos iteration scheme, applied directly to the infinite system We introduce a 
new method for the calculation of the spectral weight, that we name the Lanczos spectra 
decoding method (LSD). This new method combines the power of the Lanczos algorithm 
with some general properties of the spectral function. 

As well known, the Lanczos technique consists of a partial diagonalization of the hamil- 
tonian in the basis generated by the vectors = H"'\'$t) for n = 0, 1, A^, where |\E't) is 
a trial state. After orthogonalization of the vectors Sn, the hamiltonian reduces to tri diag- 



onal form and the resulting orthogonal basis {ei} can be iteratively calculated [jTO[- In this 
restricted basis it is straightforward to compute static and dynamical correlation functions. 
Usually, convergence is reached for N much smaller than the dimension of the Hilbert space, 
justifying the success of the method. In an infinite system this procedure can be applied 
as well, provided the Hilbert space generated by the iterative application of H to the trial 
state remains finite. This is in fact the case for the "simple" t — model if we consider 
the basis of states where the hole is located at the origin and the z-component of the spin is 
defined in any other site. The application of the hamiltonian to one of these states generates 
at most z new elements of the basis and, after A^ steps, the Hilbert space is finite having 
at most dimension . However, this exponential growth makes the problem intractable 
even for relatively small A^ . Fortunately many of the generated states appear several times 
during the expansion process of the Hilbert space which in fact turns out to be considerably 
smaller than the previous estimate. Another reduction factor can be gained by implement- 
ing the translation symmetry of the hamiltonian. Translation operators can be always used 
to move the hole at the origin of the lattice and only the corresponding basis elements must 
be stored. 

By applying the above strategy, the dimension of the Hilbert space grows much slower 



than and in fact we are able to go up to = 26 for the 2C case and = 13 for the 2D 
case with an Hilbert space dimension at most equal to ~ 12.2 x 10^. The 2C model is much 
easier to study with this method, though the properties of the spectral weight are rather 
similar to the 2D case as we will show in the following. Further details of this calculation 
will be published elsewhere P|. 



We apply the Lanczos algorithm to the t — Jz model, using as trial state \ipT > the 
exact ground state of the undoped system with a hole of definite momentum k, i.e. = 
J^ReA g^^^cr^I^) where |$) is the Neel state with spin up electrons in the A sublattice 



'L/2 

and L is the size of the lattice, where periodic boundary conditions are assumed. 
The spectral weight can be formally calculated at fixed and reads: 

A{k,u) = Im-(^'T| ^ ^I^t) 

vr u — n — 10 

N 

= Y.\{^^\^T)?5{u-E,) (2) 

i=0 

where Ei and are eigenvalues and eigenstates of H in the restricted Lanczos basis. We 
expect that for ^ the spectral function is completely incoherent in the infinite size 
limit and we know that by the NK theorem A{k,uj) is identically zero outside the interval 
ep < oj < —ep- As a result of the finiteness of our Hilbert space, at fixed A^ we get a 
sum of 5-functions in the spectral weight. This feature also appears in the exact spectral 
function of finite systems and, in that framework, an estimate of the thermodynamic limit 
is obtained by smoothing the 5-functions in Eq. with lorenzians of a given small width 
0J^: 5{uj — Ei) lm ^^^^_^g . For small 6 reasonable results can be obtained, provided 
the resolution of the energy levels becomes much smaller than 6 for large A^. In our case, 
however, we cannot reach very large A^ and a more efficient method for evaluating the 
spectral weight is necessary. The N oo limit of the spectral function may contain a 
coherent part only if in this limit some Zj = attains a finite value. We have 

verified that for both the 2C and the square lattice, all Zi tend to zero and therefore the 
spectral function is completely incoherent as expected. The spectral weight can be always 
written product of two functions: 



A{k,uj) = Z{uj)pLiuj) (3) 

where 

1 ^ 

PlM = -— (4) 
and Z{uj) is defined at discrete points: 

Z{uj) = {N + 1)\{^,\^t)\^ (5) 

for = Ei. If 74(fc, a;) is incolierent and the Lanczos eigenvalues are set in ascending order 
Eij^i > Ei, Z{uj) describes a well behaved function of u as well as the coarse grained Lanczos 
density of states (LDOS) which can be estimated from Eq. (4) as 

Ei + Ei+i 1 

- - {N + im,,-E^y 

Analogously, Z{uj) can be then interpolated linearly at energies q: 

Z(e,) = + Z,)/2. (7) 

Eqs. d^J^) are accurate to 0{1/N'^) and better fits are possible using more data for the 
interpolation procedure. A{k,uj) easily follows from Eq. (^. At the end, we can verify the 
sum rule / A{k,uj)duj = 1, as a check for the accuracy of the calculation. This is all about 
the LSD method. 

In order to test LSD we have applied it to an exactly solvable problem. If we neglect 
all closed paths in the averages of < tpT\H"'\tpT > the resulting Green's function is the 
known BR one. On the other hand, this is the exact Green's function in the Bethe lattice 
of coordination z. LSD is perfectly defined in this case. In fact, it is possible to compute 
analytically Z and Pl in Eqs- due to the simple structure of the Lanczos matrix in 

the Bethe lattice case: (ei|i7|ei+i) = \/ z — 1 for z > and (eo|if|ei) = y/z are the only 
non-zero upper diagonal matrix elements in the orthogonal Lanczos basis {ei}. 

In Fig. |I| we show a comparison of the spectral weight obtained by smearing out the 
5-functions-usual method- and by applying our scheme. We see that the agreement of 



our data with the exact results is very good and we easily resolve the first peak in A{uj) 
(dispersionless in this case) for the 2C case. By contrast the usual method cannot give 
accurate results with the available number of Lanczos steps. 

In the following, we present the results for the 2C and the 2D case, obtained using this 
new method to evaluate the spectral weight, without any approximation. We show in Fig. |^ 
the spectral weight for k = (0, 0) obtained with different numbers of Lanczos steps. For the 
2C case, our results are well converged and we find a sharp peak located at an energy close 
to the BR one Cbr = — 2tV z — 1 and, surprisingly, a second peak at energy ~ —t. In the 
2D case, the spectral weight looks similar, although the first peak is rather small. 

In ID, the exact BR solution leads only to one peak but with a divergent spectral weight 
~ at the bottom of the band. Already in the 2C case such a divergence disappears 

within the RP, as well as in our numerical scheme, which includes all closed-loop paths. In 
fact the peak in Fig. ^ does not depend much on the number of Lanczos steps. 

In Fig. ^ we see that the first peak in the spectral function has a remarkable dispersive 
feature although the bottom of the spectrum appears A;— independent. The dispersion of the 
first peak is not present neither in ID or in infinite dimension [^,|^ and the importance to 
go beyond the RP is clear even in 2D. 

The one particle density of states (DOS) can be computed either by integrating A{k,u!) 
over k or equivalently by a direct evaluation of lmG{R = 0,^;). As it is shown in Fig. 2 
for the 2C case, our results present some small deviation to the DOS in RP approximation 
especially when large closed paths are allowed. In 2D however the BR solution seems already 
quite accurate, at least away from the band tails. 

Another open problem is whether the band edge of a hole in a QAF coincides with the 
NK energy and how the spectral weight vanishes at the band edges. At finite A^ the lowest 
eigenvalue of the hamiltonian E^- restricted to the Lanczos basis is of course a variational 
upper bound to the lowest eigenstate non-orthogonal to {ipT >■ Note that the NK state is 
orthogonal to {ipr > only in the infinite system and E^o = lim Ej^j may be different from 
ep, contrary to the finite size case. Thus E^o is a definition of the one hole energy in a 
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quantum antiferromagnet and since it coincides with the smallest energy uo where A{k^uj) 
vanishes. In order to have a good estimate of this energy E^^ it is useful to know what 
are the leading corrections of the quantity Aa? = — E^o for — > oo. The way Aat 
vanishes for oo is related to the form of the LDOS at low energy. In the BR case the 

exact solution gives Pl{^) ~ e^"*^/^ (see Fig. 2). Thus, using Eqs. A^^''^ ~ ]vK^' yielding 
Aat ~ In the general case we have numerical evidence of a finite LDOS, and the same 
argument determines Aat ~ 

We have plotted in Fig.4(a) the estimated ground state energies as a function of for 
several momenta for the 2C case. Many of the estimated Lanczos energies - exact upper 
bound of the true ground state energy- are clearly below the BR energy (even for the 2D 



case not shown in the figure ||TT|). Thus, a previous suggestion that the one hole energy in a 



QAF should be close to cbr |lT2|jr3|| is not confirmedhy our numerical results. In Fig. 4(a) it 
is a remarkable property that all the extrapolated energies are very close to the NK energy, 
independent of the momentum of the hole, E^o = — 3±0.02 although the spin configuration 
is antiferromagnetic. This clearly suggests that the spectral weight is finite up to the NK 
energy. Moreover the fact that the corresponding quasiparticle weight becomes very small 
close to the NK energy (e.g. Z = 0.02 for i?26 = —2.92 and p = 0) may indicate that the 
vanishing of the spectral weight is non-analytic. Indeed a least square fit of our data for the 
band tail is more consistent with an exponential vanishing rather than with a power law. 

In order to check the estimated — energy we have also computed Ej\f for the finite 
Jz model, where due to the localization of the hole in the linear "string" potential Ejsi 
converges exponentially with and A^ = 26 gives already a very accurate estimate (within 
0.1 %) for Jz > 0.1. For smaller Jz ( 0.02 < Jz < 0.1 ) the same accuracy is obtained by 
fitting the Lanczos data with A at oc e ^ ^ -which interpolates consistently the Aat ^ 
convergence in the limit Jz = and in the finite Jz case. In the string picture the one hole 
energy for small Jz should he E = a + b^j-^ with a = cbr and 6 = 2/3 independent of z. Our 
data shown in Fig. 4(b) are clearly consistent with the BR exponent 2/3, although a = ep- 
The next leading corrections to the contribution Jz^^, are also very important. For small 



Jz in the 2C case we do not find any transition to a phase separated (PS) polaron state 
0: Eps = -3t + c(^)2/3, with c/t = 3(7r/2)2/3 = 4.05. Previous numerical work ^ is in 



quahtative agreement with our resuhs although the small size studied [jlj| was not enough 
to detect the — > smooth crossover to the NK energy, that we found in the 2C case. 

In summary, we have presented here a successful attempt to go beyond the RP for the 
hole dynamics in a QAF. A new Lanczos-type of analysis of the spectral weight enabled us 
to get very accurate results for the 2C problem and qualitatively similar ones for the 2D 
case. A clear dispersion of the main incoherent peak of A{k,uj) both for the 2C and the 
2D case was found. The variational argument about phase separation on the t — model 
(or t — J) for small Jz was based on the assumption that the one hole energy in a QAF is 



smaller than the NK energy ep by a finite amount ||T3[ . Instead our results strongly suggest 
that Eoo = ep for 0. In fact as it is shown in Fig.4(b) the one hole energy in a 

phase separated polaron is always well above the exact estimated ground state energy for 
all > 0.02. Of course numerically we cannot rule out that a phase transition to a non 
uniform phase may actually occur for unphysically small Jz- Anyway the fit of the energy 
data for small Jz extrapolate to the Jz = NK energy with a remarkable accuracy (see 
caption Fig. 4b). 

Thus the phase separation in the small J t — J model may be washed out by quantum 
fluctuations [jl5|,|l6l, especially in higher dimensions where Eps approaches ep ioi Jz Q 
with a much smaller exponent. 

We gratefully acknowledge useful discussions with E. Tosatti. This work has been par- 
tially supported by CNR under Progetto Finalizzato "Sistemi informatici e calcolo parallelo" . 
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FIGURES 

FIG. 1. Comparison of the LSD for 10, 18 and 26 (8,11 and 13) Lanczos iterations with cor- 
responding analytical results (solid lines) within the RP for the Bethe lattice with coordination 
number z = 3 (z = 4). Triangles, squares and circles correspond to the small, medium and large 
N calculation respectively. The long dashed lines and dashed-dot lines are fit of A{uj) obtained by 
the standard method (see text) with 5 = 0.05 and 6 = 0.1 respectively. 

FIG. 2. Calculated k = -spectral weight and DOS for N = 10, 18,26 (8, 11, 13) for the 2C 
(2D) model. Solid lines are cubic interpolations of the largest N, and the dashed lines are the BR 
densities of states and a guide to the eye for the 2D-^(A;,u;). The symbols for the points are as in 
Fig.l. 

FIG. 3. Calculated momentum dependent spectral function (lower part) for different k in the 
magnetic Brillouin zone and the dispersion of the first peak for the same k values (upper part). 
The bars are typical estimates of the first-peak width. For the 2C A{k, uj), the wavevector k ranges 
from from (0, 0) ( bottom) to (vr, 0) (top) with nine equally spaced values. For each k, we have 
shifted the spectral function by 0.25 successively. For the 2D A(k, iv) the fc— path in the magnetic 
Brillouin zone is shown in the inset The inset in the 2C dispersion is the free electron behavior. 

FIG. 4. (a):Plot of the lowest eigenvalues of the 2C model as a function of 1/N, the inverse 
of the Lanczos iteration number, for the same momenta shown in Fig. 3. The horizontal dashed 

2/3 

lines is the BR ground state energy, (b): ground state energies (crosses) vs Jz for the 2C model. 
The continuous line connects linearly the data and the dotted line is an extrapolation ( consistent 
with the RP 0) of the last five points for 0: E{Jz) = E{0) + a^'^^^ + bJ^, with a/t = 2.6, 

b = —2.2 and E(0) = — 2.996t. The long dashed line is the polaron energy (see text). 
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